RFoptions(seed=0, xi=0)
## seed=0: *ANY* simulation will have the random seed 0; set
## RFoptions(seed=NA) to make them all random again
## xi=0: any simulated max-staable random field has extreme value index 0
x <- seq(0, 2, 0.01)
model <- RPopitz(RMgauss(), alpha=2)
z1 <- RFsimulate(model, x)
plot(z1, type="l")
Run the code above in your browser using DataLab